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ABSTRACT 


The dynamics of a proposed marine nuclear propulsion plant, 
under the condition of sudden changes in load, are studied 
by both an analog and digital simulation. Two inherent reac- 
tor characteristics, internal reactivity feedback and the 
existance of delayed neutron emitters, which effect power 
plant stability and controllability are considered in detail. 

Comparison of the simulation results with practical data 
indicates that the Simulation represented the dynamics of a 
marine nuclear power plant with a sufficient degree of accuracy 
to be extremely useful in the study of marine nuclear power 
plants. 

Evaluation of the simulation results indicate that the 
proposed power plant is inherently stable when operating in 
the normal region. The desirability of control schemes are 
discussed and a control scheme utilizing a constant average 
coolant temperature program is implemented. This external 
control scheme significantly improved the response of the 


system. 
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Le NRO DUCT ION 


Nuclear reactor power plants for use in marine application 
differ from stationary power plants in that they must respond 
quickly to power demands which may vary throughout a wide 
range in a relatively short period of time. There are maneu- 
vering situations where the shipboard power plant may be 
required to respond to load changes of up to 50 percent ina 
few seconds. In addition to the operational differences, the 
mobile nature of the nuclear power ship coupled with the size 
and weight limitations which must be imposed upon the reactor 
shielding require that safety be a primary consideration in 
design and control of the marine power plant. 

There are several references [Harrer 1963 and Glasstone 
and M. C. Edlund 1952] that treat the nuclear reactor and 
nuclear reactor control from the standpoint of the reactor 
operating at constant output power. Schultz [Ref. 3] discusses 
the problems of reactor control from the standpoint of a 
reactor operating with varying power demands as well as when 
operating under steady power out conditions. King ~-[Ref.—4} 
considers the problems associated with nuclear power systems 
largely from the position of a mechanical engineer. 

~ this thesis investigates a marine nuclear power plant and 
is intended as an initial look into the dynamics of this type 
of propulsion system. The system dynamic equations are 


developed and the resulting fundamental balance equations are 
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Findierved by invoking reasonable first approximations. 

Initially a power plant without an external controller is 

sueiea After the initial investigation a control scheme is 
_ formulated and implemented. The plant response with controller 
is compared with the plant without controller. 

The power plant, represented in block diagram form in 
Figure (1-1), utilizes a pressurized water reactor fueled with 
low enriched uranium dioxide, heat exchanger for the conver- 


Sion of reactor power to steam and a turbime which is repre- 


sented by the changing load. The reactor and one side of the 






HOT LEG + [HR@erTie 


€OLD LEG | 
PRIMARY LOOP , SECONDARY LOOP 
| 
Fig. (1-1). BLOCK DIAGRAM OF BASIC NUCLEAR 


POWER PLANT 


heat exchanger represent the primary loop while the other side 
of the heat exchanger and turbine represent the secondary loop 


in the system. 
13 





me System, although hypothetical, is similar to the 
plant on board N. S. Savannah. Many of the initial (before 
load changes) temperatures and coefficients closely approxi- 
mate those of Savannah [{5,—-6,—-7}. In order to show the ee 
ity of the mathematical model used in this study, a qualitative 
comparison between the response of the Simulation and both 
those of Savannah and a proposed marine power plant is used 
{8}. 

The organization of this thesis follows the format oe 
described below: Chapter (2) contains a development of. the 
mathematical model used for the evaluation of system perfor- 
mance; Chapter (3) contains a general discussion of two 
inherent reactor characteristics, delayed neutrons and internal 
reactivity feedback, that affect the stability of the power 
plant; Chapter (4) presents an analysis of the analog and 
digital computer simulations of the power plant without 
external control and also contains a comparison of the results 
of the two simulations; Chapter (5) discusses some general 
aspects of nuclear power plant control and shows the specific 
control scheme used in this evaluation; Chapter (6) is a 
discussion of the concluSions arrived at by this study and 


also discusses areas of further interest. Appendices A 


through D supplement these chapters. 
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. Il. MATHEMATICAL MODEL I 


The mathematical model for the plant under consideration 
is developed by obtaining the model for each system component 
individually. The individual models together with the values 
of the specific system parameters are summarized at the end of 
this chapter. During the course of the development, it is 
necessary to make several simplifying assumptions. These 
assumptions introduce some inaccuracies in the system model. 
“However, the errors involved in these assumptions are 
Pewaely less than “the amount of uncertainty in the engineering j 
value of the coefficients used." [3] The generalized equa- 


tions given below are listed in [Ref. 9]. 


mm. REACTOR KINETICS 

The point model kinetic equations of a chain reacting 
nuclear pile have been developed in the literature [10, 11] 
and shall be used for the starting point in the development 


of the plant model to be studied. 


6 
- 7 6k-8 : 
_ (= ) n + = anes =1) 
—— oo - i. ce (2-2) 
i" Q h ee ee 


Equations (2-1) and (2-2) are state equations describing a 
nuclear reactor operating at some steady state power condition. 


The variables n, c, and 6k represent neutron density, density 
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of delayed neutron emitters and excess reactivity. The other 
terms are constant coefficients described en-page 7. 

The indexed terms in equation (2-1) and (2-2) indicate a 
multi-group representation of the delayed neutron emitters. 
This six group representation is valid for a detailed study 
of reactor dynamics, however, in the scope of this study the 
reactor is considered as a power source which is investigated 
from the terminal point of view. Therefore, without loss of 
generality, and for the convenience and simplicity (i.e., 
equipment limitation in the analog system), the kinetic 
equations are represented by a one group delayed neutron 


approximation. Equations (2-1) and (2-2) thus become: 


- 6k _ Bn _ 
n= >7-n 7 + AC (2-3) 
¢ = SB - ic (ee) 


The magnitudes of the variables n and oc are far greater 
than the magnitudes of the other system variables. In-order 
to facilitate the use of the analog computer, equation (2-3) 
and (2-4) will be normalized by dividing through by their 


steady state values n(0) and c(0). 





Sk g . 
n/nt0) = BS) 10} . re 0) a‘, = 6 a) 
— ee h 7 
sor = tf )-sur--str— a 
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Next defining: 





Cc 
N = pee and K = 2 
fok( (0) ual 


noting that in steady state: 


N=D=i1, N=D=K#=0O0; and n(0) =X c(0) 


Equations (2-5) and (2-6) thus become: 


N = (KN - N + D) (2=7) 
p= 1@ - D) (2=16) 


The term K above is the effective net reactivity existing 
within the reactor. In steady state there is a constant 
neutron population and therefore K must be zero. To establish 
K at this zero value the reactivity attributed to rod position 
must be off set by the inherent internal negative feedback 
(internal feedback is discussed in some detail in Chapter 3). 


Thus the algebraic equation for K is: 


K = K. + Ob Th + Ob Tay (2-9) 


Tp is the average fuel temperature and T is the average 


AV 


coolant temperature. The coefficients Oh and a, are coeffi- 


@iemts Of reactivity feedback. It should be noted that Op, 


and a, are negative ‘numbers. 


We 





B. HEAT TRANSFER EQUATIONS FROM FUEL AT T._, TO COOLANT 


Si 
AT Pav 


For a given change in temperature the fuel will absorb 


heat in the amount of 
de = M Ce (T. = T,) Ts aod (2-10) 


Where qe is the heat absorbed by the fuel, Mm - 1s the mass of 


Ehevruel , Ce is the specific heat of the fuel and Tt, - T, is 
the change in fuel temperature. The time rate of change of 


the heat absorbed in this manner is: 


=> ee2 _ = 
eee tr oe 82” Tg! (2-11) 
Consider one reactor fuel element of Figure (2-1). The 
temperature profile of the fuel varies along the axial length 
Z and depends upon: 
1) the heat input from the reactor. 


2) the coolant flow rate, which is constant. 


‘ 
oA 
ae 
Z Ga) ay 
5) 5p | 
{ fe 2 
tp 
Fig. (2-1). TEMPERATURE PROFILE OF 


OS BNGLE FUEL ELEMENT. 


The first of several assumptions must now be made. Assume 


that the temperature profile does not change shape or form 
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during transient conditions and that in the small signal 
analysis Ce and mp are constant. Then at any axial position 
Z the rate of heat absorption by the fuel is dependent upon 
mie time derivative of fuel temperature. Further considering 


changes from steady state where T, is a constant, it then 


follows: 
ae Semece 2 (2-12) 
The general heat transfer equation is [12] 


(2s) 


\O 
I 
~ 
<j 
= 


axial (z direction) flow can be neglected in comparison to 
radial flow. For a homogeneous right cylinder fuel element 


equation (2-13) becomes: 
k 
. a a | : 
a> = oe (x52 (2-14) 


employing the power balance for each fuel element where Power 
= time rate of change of heat 
Bowe Sago neraccd Power Power 
by the fission = absorbed + transferred 
process by the fuel to the coolant 
Considering the case where ky 1s independent of space 


variables, the (differential) power balance equation across a 


differential radius r is: 
2 


CeeTs + i 2T,i2a (2-15) 
247 Fozr 


I) 





Or rearranging 


Z 
-_ . ~P al 2 a er 
"se “me, mc (25 +22] iZe6) 








Since Ce and ky are both functions of temperature, the above 
equation is a non-linear partial differential equation. 
Following the assumption that the temperature profile does 
not Shape or form and that 2T/2r is a constant as shown in 


Figure (2-2) 





w) 
2, 
re 
Q 
Q 
KG 
—] 
© 
FUEL COOLANT 
27 _ AT _ 7 F AV 
TEMP Que iMag Lo -Iy 
2¢T 
— 
ie 
T +7 
H'-c 
r. RADIUS r TAV = 
2 1 ; 


Poa. 2-2). TEMPERATURE PROFILE ACROSS 
FUEL ELEMENT AND CLADDING. 


thus equation (2-16) becomes 


k Te c= OT 
wt ( Zo (2-17) 








The power generated is a function of the neutron population 


may. 6 hus : 
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The term | ee ee the units of sec ~L and can 
M Ce r(r, - r,) 
be related to a time constant TT): Therefore equation (2-17) 
becomes: 
a ee oe Gee (22aie 
a 7 Ty FE AV 


@. HEAT BALANCE FOR COOLANT IN THE CORE 
Under the same assumptions as stated above the rate of heat 


absorbed by the coolant as it passes through the fuel elements 


is: 


q =m_c ac (fy - FT) (2-19) 


Where qe is the rate of heat transferred to the coolant, m 
is the mass flow rate of the coolant, CA is the specific heat 


ef tae coolant, T.,, is the hot leg temperature and T is the 


C 
cold leg temperature. Employing a power balance as above 


H 


Power transferred _ Power absorbed Ht {Power santa 
to the coolant by the coclant \ eo tiheebol lex 





k T - T ie ae 
2 (So) ea, 63 a i a (2-20) 


Where Ts is the time for a unat volume of coolant to flow 
through the reactor. 
pk 


k, @ pF) 1 


° AY . 
—— en _ -2? 
Tav foe = LPeem. c 1 (Ts, TQ) a 
2 1 ca c oO 


an. 





k 


again - can be related to a time constant T... 
r(r, = ry) mo c. 2 


Thus: 
° cael: 
i 
2 
Rearranging 
Fe 
AV To 
where 
ats 
'_ = 
Le T + 
Oo 


D. HEAT TRANSPORT 


There are time 


(2222) 


— + —— (2-23) 


DELAYS 


delays associated with the transfer of heat 


from one portion of the plant to another. In developing the 


model for the time delays the following assumptions must be 


made: 


1) there is no heat loss in the piping between the reactor 


and the boiler. 


2) there is no mixing of fluids in the system piping. 


1. Time Delay From Reactor to Boiler 


The boiler inlet temperature, T 


of the reactor outlet temperature T 


BI! is a delayed version 


H°’ 


Eo (tT) (2-24) 


BI H 


Lia 





Where Tar is the boiler inlet temperature and T3 is the 


transport time between reactor and boiler. Rearraning 


Ti = Tor (t + T3) (2-25) 


Expanding in a Taylor Series 


Le = Ter - T T By + Swmees: (2-26) 
Selving for Tae 
post @i t+2aPp 47) (2-27) 
BI T3 BI AV ‘@ 


2. Time Delay From Boiler to Reactor 


Following the above reasoning 


te —s—— (2 - 29T_. + TT...) (229) 


Where T. is the average temperature in the boiler. 


E. HEAT BALANCE FOR HEAT EXCHANGER (BOILER) 

The development of the model of a counter flow heat 
exchanger in which there is a change of phase (water to steam) 
in one flow path is a difficult task. It is not the intention 
of this chapter to present a rigourous derivation of the equa- 
tions used in the system model, but rather it is to give 
insight into the physical relationships from which these 
equations were obtained. Therefore a more heuristic argument 


will be used in the development of the steam generator model. 
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TEMPE RATURE 





DISTANCE 


Fig. (2-3). TEMPERATURE PROFILE FOR COUNTER 
FLOW HEAT EXCHANGER. 


There are several additional assumptions which must now be 
made: 
1) There is no heat lost or gained with the surroundings. 
2) The specific heat of the individual fluids remains 
Cons cane. 
3) The overall coefficient of heat transfer (k,) is 
uniform through the heat exchanges. 
4) From Figure (2-3) the mean temperature difference 
throughout the heat exchanges Th =e as i 
i  Pesinany-—Sude 
Consider a small section of heat exchanger having a 


Eaonsteresasea GA. Then the incremental heat transfer is: 


= = —-2 
dQ kK, aA at = Mm aT. (2-29) 


integrating yields 


a= kK, ALT, = Te] — in edb - TJ (2-30) 





Since ky, Ee, Cus and m, are all constant, equation (2-30) 
reduces to 


2(T, = Ts) = E(T3y = T 30) (231) 


Where To is the steam temperature and A is the area of the 
heat exchanger. 
2. Secondary-Side 
Energy balance in the secondary-side consist of the 
amount of power transferred to the secondary, KjA(T, = To), 
equal to the power absorbed within the boiler plus the power 
(eaieh ws transferred to the turbine. To bring this relation 
into a linear, ordinary differential equation some additional 
assumptions are required. These are: 
1) The thermal capacity of the boiler metal, steam 
and coolant can be lumped. 
2) The power demand on the boiler can be measured as 
a fraction of normal load. w will represent this 
condition with ~ = 1 representing the normal load. 
3) The heat exchanged in the boiler occurs at a single 
point within the heat exchanger. 
Therefore the flow of heat delivered to the turbine may be 
represented as the product of the throttle opening, represented 


by & and the steam pressure. Since the steam pressure 1s 


i ied 
Ze 

upon the steam temperature, the term Sa T. can be used 
S 

to represent steam pressure. Thus the heat balance for the 


secondary side becomes: 


e S i 
K,A(T, - Tg) = (yc +m. c.) T. + Kp WV op T (2-32) 





Where M and Ce are the mass flow rate of the steam and the 


pecerric heat of the stGam respectfully. Rearranging 


kiA Kea) 2P 
ie ~ in : + in) © [T, ~ Tg] u m = +m_c_) a the 
Cc So s ss es Cc s ss S 
(2-33) 


; 1 2P. 
Te =~ § (Th + Ts + KY og} aa) 
3) S 
a ee AL MODEL WITH SPECIFIC SYSTEM PARAMETERS 
The development of the specific equations for the system 
is carried out in detail in Appendix A, they are Summarized 


here for convenience. 


1. Reactor Kinetics 


N = 6.4 (KN-N+D) (2-35) 
D= .1 (N-D) (ere) 
ake Me otexgt0 i. - 3.14 x lol“ ota, (3-37) 


where N(0)=1, D(0)=1 and K,(0)=24.7275. 


2. Heat Transfer From Fuel to Coolant 
(T,. ~ T ay) + 30.0 N (2-38) 


where Tp (0)=565°F. 


3. Heat Balance for Coolant 


Tay = 10 ts = 50 Tay + 40 ie (2-39) 


— °o 
where Tay (0) =505 F. 
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4. Time Delays 

(T ee et + T ) (2-40) 
me fe] 

where T 37 (0) #520 jee 


esas _ ¥ 
To - 2 (T Tar/? 


wo 
AB, 
Lan neal 


(2-41) 


where T (0) =490°F. 


2. Heat Balance for Steam Generator Secondary-Side 


= 0667 + Ol0G6s v) iB (2-42) 


where T, (0) =565°F. 


G. HETHOD OF SOLUTION 

Because of the multiplication of the time dependent 
variables K and N in equation (2-7) and because K is a function 
of N the entire set of dependent equations which describe the 
plant dynamics is non-linear. Since (in general) non-linear 
differential equations do not lend themselves readily to 
analytical solutions, computer simulation will be used to 


solve these equations. 
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III. INHERENT REACTOR CHARACTERISTICS 


There are several important characteristics of nuclear 
power plants which greatly effect the stability and control- 
woerity Of the System. “The two most signii£icant factors 
relating to the system performance are the existance of 
delayed neutron emmiters and the negative internal reactivity 
feedback associated with temperature changes within the 
system. Both of these phenomena are discussed below in some 


@etail. 


A. PROMPT AND DELAYED NEUTRONS 

At the instant of fission there are neutrons (prompt) 
which are produced instantaneously. Fortunately, for control 
purposes, a small portion (8) of all neutrons produced during 
the fission process are delayed in time [l, 2]. These delayed 
neutrons are the result of the decay of fission products. As 
noted previously these delayed neutrons are grouped in six 
distinct groups with different delay constants (i). Analysis 
of the effects which delayed neutrons exhibit on the reactor 
kinetics can best be shown by the evaluation of the kinetics 
equations (2-3) and (2-4) with and without the delayed 
neutrons considered. 

The analysis of the prompt and delayed neutron response 
which follows is a study of the reactor operating at a power 
level low enough so that feedback effects (discussed later 


in this chapter) are negligible. An analagous situation in 
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a power plant would be in the start-up power range of the 
memccor. 

In this evaluation the input to the reactor is in the 
form of a perturbation to the steady state reactivity k. The 
input was limited to a maximum value (8) of (0.0064), for 
when 6k exceeds 8 the reactor is in a response condition 
described as "prompt critical." "Prompt critical" refers to 
the fact that the nuclear fission chain reaction can be main- 
tained by means of prompt neutrons alone. If this condition 
occurs neutron density, and hence power increase rapidly 
from the instant of input, making the reactor difficult to 
control and may cause what 1s referred to as a start up 
accident. The "prompt critcial" operation is avoided in 
peactice. ° 

1. Kinetic Equations Assuming No Delayed Neutrons 

Under the assumption of no delayed neutrons, equations 

(2-3) ‘and (2-4), which describe the reactor kinetics in steady 
state with a one group delayed neutron approximation, reduce 


to the single equation: 


ne — (sah 


Soigeng for n 


bk y 


ee nn eee (3-2) 


Dividing by n(0), substituting the system parameters found in 


Appendix A, and assuming 6k = 0.003 equation (3-2) now becomes: 


eee (3-3) 
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Zee aecor Kimeele Equa ons with 4a One Group Delayed 
Neutron Emitter 


= Under the one group delayed emitter assumption, 
equations (2-3) and (2-4) are unchanged. With the substitution 
of system parameters from Appendix A and from the solution of 
these equations carried out in Ref. [11], the solution of 


these equations for N is: 


+0.0884 t 3.4 t 


N= 1.88 e - 0.88 e | (3-4) 


3. Comparison of Results 
The effect of the delayed neutrons is readily seen 

by the comparison of equations (3-3) and (3-4). With no 
delayed neutrons the reactor neutron population has increased 
by a factor of 8,100 in 3 seconds, however when delayed 
neutrons are present the reactor neutron population has only 
increased by a factor of 2.5 in 3 seconds. Plots of N versus 
time for both of these situations are shown in Figures (3-1) 
and (3-2). As time passes (1i.e., time greater than 10 seconds) 
the first term of equation (3-4) becomes dominant and the 
response will follow a positive exponential increase. However, 
in the time span of interest (a small time after the distur- 
bance) the affect of the delayed neutrons is to greatly reduce 
the rate of change of neutron flux. It is the natural occur- 
ance of these delayed neutrons which makes it possible to 
operate a nuclear reactor with the degree of safety required. 
Even as small a fraction as 0.64 percent of the total neutron 
population acting as delayed neutrons has the effect of making 


the entire reactor control problem greatly simplified [2, 3]. 
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It is considered pertinent to re-emphasize that the 
above .discussion of prompt and delayed neutrons was conducted 
for a reactor in the startup range. The response of a 
reactor operating in the normal power range (where reactivity 
Peedback is strong) is significantly different than that of 
the reactor operating at low power. This difference will be 
Gemonstrated later in this chapter. 

’ It should be pointed out however, that the delayed 
neutron phenomenon also has an undesirable facet. While it is 
readily obvious that it restricts the growth of the neutron 
population for positive 6k inputs, it also restrains the 
reduction rate of neutron population which would result from 
a negative 6k input as the power demand is reduced during 


normal operation or during shut down. 


B. INTERNAL REACTIVITY FEEDBACK 

In the operation of a nuclear reactor the power generated 
is directly proportional to the neutron population [1]. In 
order for the neutron population to remain constant at some 
steady state power level, equation (2-7) must be equal to 
zero. For this situation to occur, K, the net internal 
reactivity which is the sum of both the reactivity attributed 
to the control rod position and the reactivity due to internal 
feedback, must be zero. Changes in the feedback reactivity 
depend upon temperature variations and fission-product--— 
poisoning which are discussed below 

Gogmiemre the stability of the reactor the overall coeffi- 


cient of reactivity feedback must be negative [13]. 


38 





1. Temperature Effects 


When there 1s a power change in a nuclear reactor 
there is considerable additional energy released. Since it 
is not possible to transfer all of the liberated energy to 
the load, there are temperature fluctuations throughout the 
reactor and associated coolant. 
Temperature variations affect reactivity in two 
distant ways: first, the mean energy of the thermal neutrons 
and hence their nuclear cross-section vary with temperature; 
and, second, the mean free path length and leakage probabilities 
are functions of temperature [2, 3]. These effects result in 
a net negative reactivity. A reactor power plant of the type 
described above, having a negative temperature coefficient of 
Geaectivity will tend to be self-controlling in response to 
load changes and limited external reactivity variations [14]. 
In hetrogeneous reactors (the type under study here) 
the fuel and the coolant are at different temperatures and at 
different physical states. Therefore, there will be specific 
effects associated with both the fuel and coolant [15]. 
Temperature effects upon reactivity manifest themselves 


in temperature coefficients of reactivity. (a. and a). 


EF 
Ze Fsseon-Product Poisoning 


When a reactor is operating at power certain fission 
products, notably Xenon 135 and Samarium 149, are accumulated. 
These products have large nuclear cross sections and absorb 
the excess neutrons, which are vital in maintaining the chain 


reaction. The rate of formation of these products is 
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dependent upon the rate of fission occuring at the specific 
Bewer level. Eventually the rate of formation and the loss 
of the absorbing nuclei become equal and equilibrium exists. 
However, when the reactor is shut down, these substances, 
which result mainly from the decay of other fission products, 
continue to grow and may reach their highest concentration 
several hours after shut down occurs. Since the production 
of these fission products continues after shut down, it is 
entirely possible that there may be enough negative reactivity 
remaining to restrict the power build up of the reactor until 
these products decay to a level where the oe positive 
reactivity in the control rods can overcome their adverse 
reactivity effects. 

The above fission-product phenomenon does not essent- 
ljally affect the performance of a reactor operating at power. 
For this reason these effects are not considered in the scope 
@e this study. 

Se ceeinwepnalsReactiviity 
The net reactivity in dollar units, due to both rod 


position and to temperature can be expressed as: 


C. TOTAL EFFECT OF DELAYED NEUTRONS AND NEGATIVE REACTIVITY 
FEEDBACK 


The total effect of the delayed neutrons and negative 
Peaeelvity Can best be illustrated in Figure (3-3). These 


plots are for the specific plant under study. The power plant 
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disturbance is a step input of magnitude K = 0.469 dollars 
fmee., Ok = 0.003). 

The input causes the reactor's neutron density to rise 
sharply from a steady state unit value to a value of 1.6. If 
there were no delayed neutrons this neutron density would 
rise faster and obtain a larger value. The negative reactivity 
feedback acts rapidly and limits the excursion of the net 
reactivity to a maximum value of 0.38. After approximately 
60 seconds, the reactivity feedback has adjusted the internal 
reactivity to zero. The steady state effect of the disturbance 


is to raise the neutron density to a level of 1.02. 
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IV. DIGITAL AND ANALOG SIMULATION OF POWER PLANT 
WITHOUT EXTERNAL CONTROL 
A. DIGITAL SIMULATION 

The system equations were programmed and solved on the IBM 
360-67, using the Naval Postgraduate School Library Program 
(DRKGS) for the integration. This integration method utilizes 
a forth-order "Runge-Kutta" solution with double precision 
arithmetic. 

Time response of the system to load changes of +40 percent 
are shown in Appendix B and are tabulated for all experimental 
load changes at the end of this chapter. 

1. Load Changes 

The system was subjected to standardized step load 
@eamees distusbances of + 20 percent and + 40 percent. These 
load changes were selected because a 20 percent load change 
represents a change which a maritime power plant should be 
able to respond to quickly, and a 40 percent load change 
represents a feasible change which could be encountered during 
an accidental situation. 

2. Results | 

The results of the simulation of this hypothetical 
system were compared against the time response of marine power 
plant designs of similar size [6, 8]. It was found that the 
Simulation used in this study compared well with the referenced 


systems. The steady state (after the load distrubance) values 
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of the system variables coincided well with their counterparts. 
The simulated system demonstrated a slightly more oscillatory 
Pesevior @imem the actual and propesed systems referenced. 
This oscillatory behavior is attributed to the simplified 
system model used in this study. Even with Biagerieaigl oscil- 
latory differences, the simulation represents the actual 
metem with a sufficient degree of netbaet that Lt 1s ‘quate 
Suitable for an initial study of marine reactor power plants. 
3. Discussion of Results 

a) The power plant, without benefit of an external 
control system, was shown to be stable to all subjected load 
changes. The system displayed a second order well damped 
response to load changes. ‘ 

b) Temperature changes, including overshoots, are not 
severe and should not cause thermal stress problems. 


c) System response, although sluggish, is fast enough 


for use aS a maritime power plant. 


B. ANALOG SOLUTION 

The analog simulation of the system was carried out on 
Ene Comeor CI-3000 computer. This computer has a + 100 volt 
dynamic range. The CI-5000 computer is interfaced with the 
XDS-9300 digital computer. This hybrid pair will be used 
later in the control section of this study. 

1. Prompt-jump Approximation 

Figure (B-l1) of Appendix B shows that the time rate 

of change of N is not large. For ease of analog simulation 


(i.e., to reduce the number of high gain amplifiers required 
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in the simulation) equation (2-7) can be modified as follows: 


a. 


a KN -N+D (4-1) 


making the approximation that N = 0 and solving for N 


This approximation is called the prompt-jump approximation 
Since any change in K is reflected instantaneously in N. The 
system equations using this assumption were simulated on the 
IBM 360-67 and it was found that this approximation does not 
introduce errors greater than 2 percent in the system response. 
2. ssVoremevantables to Computer Variables 

The system variables such as temperature and relative 
neutron population must be related to analog computer variables 
which are voltages. The system equations and initial condi- 
tions in terms of computer variables are developed in Appendix 
C. These equations are summarized below for convenience: 


a. Reactor Kinetics 





- 5D _ 

ee a (4-3) 
SS Oil ese (4-4) 
K =o oko —- 08471 T. — 0.949°T (4-5) 
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where 


N(O) = 10 volts 
D(O) = 10 volts 
K, (0) = 35.36 volts 


Heat Transfer From Fuel to Coolant 


Loe oes i =f) + 0.3 N 
Where 
T, (0) = 94.17 volts 


T ay (0) — "G4 .17 Vourts 


Time Delays 


Tar = -0.5 [Th - 2 PAV + TOI 
me -- _- i alesse 
T. = 0.5 [ To a z To “ x Tay! 
Where 

Thy 00) = 86.67 volts 

T 0) ="ol-07. VOLES 
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(4-6) 


(7) 


(4-8) 


(4-9) 


(4-10) 





e. Heat Balance for Steam Generator Secondary-Side 


T. = [0.133 Tar — O,aez4 Ty = ORO OZ" aa) To] (4=11) 


ET =3979H117 volte 


The analog computer diagram for the solution of these equations 
is shown in Figure (4-1). 
3. Comparison of Analog and Digital Solution 

Below is shown two tables comparing the results of the 
digital and analog computer solutions. Because of the inherent 
resolution quality characteristics of the digital computer, 
the digital solution will be used as the base for determining 
errors. Settling time is defined as the amount of time 
required for the magnitude of K(t), the net internal reactivity, 
Peeeeeetunero a Valuc of 4 0.005 after the applied step load 
disturbances. 

From the tabulated ee ee 1t 1S apparent that the 
analog and digital solutions are in such close agreement that 
the analog solution is adaquate for the remainder of this study. 
Henceforth the digital computer is used as an auxiliary compu- 
tational tool to solve specific problems but is not used in 
the evaluation of the overall simulation. 

The value of the system variables obtained from the 
analog simulation were obtained through Analog to Digital 


conversion and printed out from the XDS-9300 digital computer. 
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Fig. (4-1). ANALOG COMPUTER DIAGRAM FOR BASIC POWER 
PLANT. 
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TABLE (4-1). COMPARISON OF ANALOG AND DIGITAL SOLUTION 
BOR A SIEP LOA entice Or yz0e PERCENT. 
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TABLE (4-2). COMPARISON OF ANALOG AND DIGITAL SOLUTION FOR 
A STEP LOAD CHANGE OF 40 PERCENT. 
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V. CONTROL 


The purpose of an external control scheme for a nuclear 
power plant must be to set up the desired steady state 
Operating conditions, and to restrict transients in vital 
plant variables to be consistant with design limitations. 

The above purposes should be obtained without any significant 
reduction in the inherent stability characteristics of the 
plant without external control. 

Examination of the basic plant transient responses in 
Appendix B indicates that the system under study is very 
stable when subjected to load changes. It is however, 
desirable, as pointed out Seow to limit the excursion of 
certain plant variables from their steady state normal value. 
In a plant with a large degree of inherent stability direct 
control of plant variables may be attempted. If a plant has 
poor inherent stability, neutron-level control should also be 


used [3]. 


A. STEADY STATE PROGRAMMING 


The pattern that the plant temperatures, pressures, and 
flow rates assume as a function of power-demand is oftened 
referred to as a program [4]. A specific control strategy 
can be devised to obtain the desired program. The two most 
Significant programs for discussion are the constant Tay and 


senstant To programs. 


/ 
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b. Censtant T 


AV 
e ie 
TEMPERATURE T 
OR AV 
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POWER DEMAND 


Bag. (5-1). SKETCH OF TEMPERATURE AND PRESSURE 


PROFILES FOR A CONSTANT Tav PROGRAM. 


Figure (5-1) is a graphical representation of the 
constant Tay program. At zero power demand, all significant 
plant temperatures (except T.) are at the same value. As the 


power demand increases, the steam temperature and hence 


pressure decrease, while Ty and To vary in such a manner as 


to keep T constant. 


AV 
2. Constant T. a 

a 
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fo-2) SKETCH OF TEMPERATURE “AND PRESSURE 


Phi. 
PROFILES FOR A CONSTANT To PROGRAM. 
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Again at no load all pertinent temperatures are equal. 
As the load increases the temperature and the pressures follow 
the pattern as shown in Figure (5-2). 


3. Comparison of Constant Tay and Constant To Programs 


The constant Tay Program is the one preferred by the 


reactor. The principle advantage of the constant T program 


AV 
for this type reactor is that the primary coolant pressure 
also remains constant. Therefore the need for an external 
pressurizer in the primary circuit is reduced. 

One additional advantage occurs if the plant under 
study had been fuel with an uranium alloy vice an oxide of 
uranium. For a metal alloy fuel the feedback of internal 
reactivity associated with changes in fuel temperature is 
insignificant when compared to the feedback associated with 


changes in average coolant and therefore can be ignored [3]. 


Bor this case equation (2-9) becomes: 


e—- — + co Tf (5-1) 


Since in stead stake K must equal zero, it is apparent 
from equation (5-1) that there is no need for control rod 
motion to maintain Tay constant. However the plant under study 
is fueled with uo, and changes in fuel temperature have a 
strong affect on net reactivity. Examination of equation 
(2-9) now indicates that there 1S a requirement to change 
control rod position to keep the average temperature constant. 


The advantages of the constant T. Program are apparent, 


S 


faem Figure (5-2), in the design of the secondaxy loop of the 
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plant. A constant To program also provides a constant steam 


pressure which allows for an optimum design of the steam 
plant with conventional controlling devices. 
The two primary disadvantages of the constant Te 


program are first, allowing T to vary over a large range 


AV 


(required to keep T, constant) may cause serious problems in 


> 
the primary coolant pressure system and second, when the plant 
is changing from one steady state power level to another, the 
plant must fight the tendency of the coolant negative temper- 
ature coefficient to hold Tay constant. 

With the exception of Nautilus, which uses a constant 
To program, all Naval plants use a constant Tay preogzam [16 9" 
The N. S. Savannah utilizes a control plan for a constant Tay 


program with a plus or minus 3°F dead zone about the normal 


load Tay value [17]. 


Be WEoLGN OF THE EXTERNAL CONTROL SYSTEM 


Consistant with the advantages of the constant T y Program, 


A 
the simulated power plant will be controlled under a constant 
Tay Program. As noted in section (5A-3) control rod movement 
is required in order to obtain the desired program. The next 
decision to be made is at what rate may reactivity be externally 
introduced into the reactor. The answer to this and other 
design questions are discussed below. 
Mroggeectivity Rate 
The one basic underlying consideration which dictates 


the design of a control system for a nuclear power plant is 


safety. As shown in Figure (3-3) the power plant operating 
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at normal power can withstand a step reactivity input of 
magnitude OK = .469$. This however is not the case for a 
acon operating at reduced power levels where the negative 
reactivity feedback effect is reduced. Therefore in order to 
provide a maximum safety margin to avoid startup or low power 
accidents input reactivity rates will be restricted to a value 
which can be handled by a reactor See oG at lower power. 
The solution to equations (2-7) and (2-8) for ramp input of 
reactivity, K(t), does not give a valid solution to the 
problem because there is no feedback considered. Externally 
introduced reactivity rates are determined largely from 
empirical data on the individual plants involved. A reasonable 
reactivity rate for a plant of this size is 10¢/sec (0.1 S$/sec) 
i). 

It 1s feasible that this rate of reactivity may be 
controlled by a manual means. However, it was desired that 
the system be fully automated therefore an automatic controller 
was used. 

2, ®@bomatic Controller 

The automatic controller consists of a hybrid combina- 

tion of an analog control rod positioning system and a digital 


Soieem Lox comparison of T plus the time derivative of T 


AV AV 
with the reference value of 505°F. It was decided to use the 
+ 3°F dead band error detector of N. S. Savannah as a valid 
method to obtain a discontinuous controller. The principle 
advantage of a discontinuous controller is that it utilizes 


the reactor inherent stability characteristics to reduce wear 


on the control rod mechanisms. The dead zone, although it 
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Mmst be consistent with the design temperatwre (i1.e., pressure) 
MeuctUatlens in the primary loop, discrimima@ies against noise 
on the Tay Signal. A block diagram description of the entire 
Simulation, plant plus controller, is represented in Figure 
(5-3). The design of the control rod positioning system and 
a discussion of the hybrid simulation along with transient 
response curves for a step load input of +40 percent of the 
system with external control are shown in Appendix D. 
3. Evaluation of External Controller 

The external controller was implemented and incorpor- 
ated into the simulation. The system was subjected to step 
load disturbances of plus and minus 20 and 40 percent. The 
improvement in system performance can be seen from tables (5-1) 


mai 5-4). The reduction of the excursions of T both 


AV’ 
transient and steady state, and the reduction of settling time 
were considered significant. 

The plant was subjected to step load changes of up to 
80 percent (an unrealistic change). It was found that the 
system was stable and that the system response improved up 
to step load inputs of 60 percent. Above 60 percent, the 
system response deteriorated and eventually went into a limit 
cycle when the step load disturbance became greater than 75 
percent. The deteriorated performance and eventual instability 
of the system were attributed to the fact that the derivative 
Of T Vy became very large and began to work in the opposite 


A 


direction of the Tay excursions rather than in normal anti- 


cipatory mode of a position plus derivative feedback controller. 
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TABLE (5-1). COMPARISON OF SIMULATED POWER PLANT 
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iMerconteroller was constdered adequate for all 


reasonable power changes. 
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* VI. CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions and recommendations are offered 


as a result of the investigation contained in this report: 


A. CONCLUSIONS 

1. <A properly designed pressurized water nuclear reactor 
(1.e€., one with a negative temperature sOMEticient OF 
reactivity) possesses a great deal of inherent stability when 
the associated power plant is operating in the normal power 
range. It was observed that the negative temperature coeffi- 
cient of reactivity served to stabilize the reactor for dis- 
turbances in both load and input reactivity. 

2. The naturally occuring phenomenon of delayed neutrons 
greatly simplifies the control problems associated with 
nuclear reactors. The effect of this phenomenon is more 
evident at the lower power levels before the temperature 
reactivity feedback effect has become affective. 

3. The control of a stable nuclear power plant operating 


in the power range under a constant T y Program with constant 


A 
coolant flow rate is a relatively simple task. If sufficient 
Stability exists in the reactor then direct control of a plant 
parameter (Tay) can be attempted. 

4. If step load changes greater than 60 percent, for which 
the above control system is not suitable, are considered 


feasible a comparator could be installed in the differentiating 


circuit to disregard the derivative if it becomes too large. 


56 





B. RECOMMENDATIONS 

Further workers in this area should comsider the following 
@acets of reactor systems. 

1. The above study was concerned with the reactor power 
plant dynamics when excited by both load disturbances and 
changes in reactivity investigated separately. Another area 
of interest is the response of the reactor under the condition 
Of load changes and externally induced rea@etivity changes 
(both aiding and opposing) occuring simultaneously. This 
Situation might represent an accidental condition. 


2. The Tay PLegram used EO control the) power plant does 
Ty + T, 
not consider the maximum value to Thy Since Tay = rr 


the value of T,, might become too high for design considera- 


tions. To institute a control program with a constant Tay 
and also limiting T,, to some maximum value would make the 
requirement of a variable coolant flow rate necessary. This 


added requirement would necessitate a new mathematical model. 
Under the assumption of a constant flow rate, the time coolant 
spends in the reactor (Ty) is a constant. However, in the 


model for a plant where a maximum value of T,, must be 


H 


considered To would be a function of Th 

3. More sophisticated control problems, than reactor 
control at operating conditions, exist in the nuclear power 
plant control field. These are the problems associated with 
the reactor in the startup range and problems involved in the 


mechanical control of the various plant components of the 


secondary loop. The problems in the control of the secondary 


Si 





loop result from the fact that when the reactor is controlled 
under a constant Tay program, the steam temperature and 
pressure vary Significantly over the operating range of the 
plant. It is recommended that further investigators consider 


the control problems of the secondary loop. 
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APPENDIX A 


SPECIES SYSTEM EQUATIONS 


In the development of the generalized system equation, 
equations (2-28) and (2-34) are expressed in terms of the 
variable T,- T is not a state variable of the system. The 
equation for the primary side of the heat exchanger (251) Can 
be solved for T. and this value substiuted into the equations 
(2-28) and (2-34) to place them in a more suitable form. 


This yields: 


= a ec a E-l op 
Cc ‘ C e+] E+l “BI 


The system parameters for the hypothetical plant are stated 
below. It should be emphasized that these parameters do not 
represent any actual real system but that the parameters used 
here are obtained from various sources in the literature [5, 
6, 7, 8, 9] and represent typical values which might be 
associated with a plant the size of the power plant of N. S. 
Savannah. 

1) Fraction of delayed neutrons B = 0.0064 

2) Average delay constant of the 6 group delayed 


neutrons KX = 0.1 
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3) 


4) 


>) 


6) 


7) 


8) 


9:) 


Mean neutron life time 2 = 10” sec 


Steady state temperature 


A) Average fuel temperature T,(0) = 565°F 
B) Average coolant temperature T ay 69) = 505°F 
C) Boiler inlet temperature Thy 69) =e 2 oe 
D) Steam temperature T, (0) = 475°F 
E) Cold leg temperature T. (0) = 490°F 


Temperature coefficient of reactivity for the fuel 


(UO.) 


tees -beS7T x ion $/°F 


Temperature coefficient of reactivity for the coolant 


(pressurized water) 


a. = -3.14 x 10° $/°F 


Time constants for the various heat transfer equations 


A) tT) = 2.0 sec. 
B) ee 0.1 sec. 
C) T3 = 2 0usee. 
D) ne 2.0 Sec; * 
E) ic. > Se OeSeens 


In the region of interest the steam pressure-steam 
temperature relation is approximately constant 


2P. 
sa = 5 psi/°F 
Us 


The values of K, (0), Ayr K tT, were determined 


red ATs 


by setting the equation in which these terms appear to 
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zero (steady state) and solving for the unknown 


Varmncables. The results are: 


A) K, (0) = 24.7275 §$ 

B) Ay = 30.0 

Cymer =O 

Dae el 1saaye ton 

E) To = 0.05 sec. 

F) To = 0.02 sec 

*T 31 Tar and Te were obtained by engineering approxi- 


mation and are substantially less than those mentioned 


for a stationary power plant [9]. 
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APPENDIX B 


TIME RESPONSE OF POWER PLANT 


As a representative example of system time response 
Figures (B-1), (B-la), (B-2) and (B-2a) of this appendix show 
the time response of selected system variables to a step load 
change of +40 percent normal load. Figures (B-1) and (B-1la) 
are solutions using the digital simulation. Figures (B-2) 
and (B-2a) are solutions utilizing the analog simulation. 
Pertinent values of the variables are shown adjacent to their 
time response. The system response for other experimental 


load changes are similar in nature to those shown. 


62 





1.0 


l-o 


K(MAX) = .119 §$ 


40 80 120 160 secs 
NEES REACTIVITY VERSUS TIME 


N (MAX) = 1.38 
N(SS) = 1.37 
40 80 mo 160° secs 


NEUTRON POPULATION VERSUS TIME 


WwW = 1.04+0.4u (t-20) 


40 80 120 160 secs 


LOAD CHANGE 


Fig. (B-1). TIME RESPONSE OF SIMULATED POWER 
PLANT TO A 40 PERCENT STEP LOAD CHANGE 
(DIGITAL) . 
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T_ (MIN) = 455.3°F 
YS °F s | 
\ T.(SS) = 456.4°F 
0 40 80 120 160 secs 
STEAM TEMPERATURE VERSUS TIME 
—. ae | 
“2 + T.(MAX) = 580.1°F 
T.(SS) = 579.9°F 
0 40 80. Lae 160 secs 
FUEL TEMPERATURE VERSUS TIME 
JOS" 
| = ° 
\ Pay (MIN) 496.3°F 
po O 
P.y(SS) = 497.5°R 
0 40 80 120 £160 secs 


AVERAGE COOLANT TEMPERATURE VERSUS TIME 


Big. (B-la). TIME RESPONSE OF SEMUBATED POWER PEAT 
TO A 40 PERCENT STEP LOAD CHANGE (DIGITAL). 
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LOAD CHANGE 


Fig. (B-2). TIME RESPONSE OF SIMULATED POWER PLANT 
TO A 40 PERCENT STEP LOAD CHANGE (ANALOG). 
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= ° 

fee eae eee oa), To) = 450.0°F 
ae | eee ee eee 

ee a ee ee] sot.(ss) = 457.0°R 
Be Jee See Se 
yosor LL ER ee 

T (MAX) = 580.4°F 

ieSUSIS) Sk) cals 

— ° 

Tay (MIN) = 495-0 24 

— ce) 

Tay (SS) = 949 720 2 





AVERAGE COOLANT TEMPERATURE VERSUS TIME 


Bag. (B-2a). TIME RESPONSE OF SiMe P FOveR Piald 
TO A 40 PERCENT STEP LOAD CHANGE (ANALOG). 
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APPENDIX C 
PREPARATION OF SYSTEM EQUATIONS 
FOR USE ON ANALOG COMPUTER 

Equations (222) Mand “(2=36) through (2-42) describe the 
power plant in terms of system variables. These system vari- 
ables have units such as degrees Fahrenheit and dollar units 
of reactivity which must be related to analog computer vari- 
ables which are voltages. In addition to the conversion of 
units, the analog voltages have a maximum value of 100 volts 
SO Scaling of the variables may also be required. 

To accomplish this conversion a scaling factor coefficient 
me introduced. The relation between computer variables and 


system variables then becomes: 


Tp = xX T 


and so forth. 
The individual x's are equal to the maximum expected value 
of the system variable divided by the maximum desired voltage 


for the computer variables. For example: 


-~s 
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assuming T,, (MAX) = 600°F and Ti (MAX) = GIG volts 
oi 
‘ Xmp = cl = 6°F/volt 


A summary of the scaling coefficients is as follows: 


xy = fF yore Ps 

Xp = vee — 

x, = Zoo Lt 
Xep = 7 S7vomt 

<w0e—.] vort = 

cee ts “Pay “TC “men 7 > oat 


The system equations are now altered by replacing the system 
variables with their corresponding computer variables and 
conversion factors. The initial condition of the differential 
equations are determined as follows: 


T,, (0) 





JN 


In a final form the equations now become: 


1. Reactor Kinetics 


Al 
it 
WW 
nn 
~ 
1 
ean) 
= 
~I 
jy 
KH 
] 
es) 
\O 
=> 
NO 
Z 
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2. 


oe 


4. 


De 


N(O) = 10 volts 
P40) = 10 vodses 
BO) = 35.36 volts 


Heat Transfer from Fuel to Coolant 


tae — Od (T., = Desi ta 


F Tay 
Where (0) = 94.17 volts 


Heat Balance for Coolant 


—S ES jee SSS EE 


Tav = -50 Tay iT Tp + 40 To 


Where T 


ay £9) = 84.17 


Time Delays 


T =~ Orso (-Ths + oT ay = Ta) 


Where TR 7 9) =" 86.607. Vola 


TA (0) = 81.67 volts 


Heat Balance for Steam Generator Secondary-Side 


i = (OR133°rR, — 132 SsorvC ema 
Where T. (0) = 79,17 volts 
wo = 5 volts 
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APPENDIX D 


CONTROLLER SIMULATION 


Analytical design of the external controller could have 
been approached in one of two ways. First a small signal 
linearization of the plant dynamics could have been attempted. 
This method was not valid since two of the plant variables, 

N and D, change in almost direct proportiom to load changes. 
For example a 40 percent load change causes a 40. percent 
change in N and D. This change in magnitude of the variables 
was outside the range of small signal application. Second a 
describing function could have been obtained utilizing sinu- 
sodial inputs [18]. The describing function technique is a 
laborious method of obtaining insight into the absolute 
Stability limits of the system but gives little insight INS) 
the controller design problem. In this initial study the 
controller design and evaluation was carried out through 
Simulation. If a more detailed study had been undertaken, a 
describing function would have been necessary to verify the 
Simulation results. 

The simulation of the external controller was accomplished 
Py Wise Of a trial and error technague utilizing the hybrid 
Gemputing feature of the XDS-9300 and Comeer CI350006cemputers. 

The analog voltage representing Tay was sampled at a 
frequency of 10 samples per second (this sampling frequency is 


well in excess of that required by the sampling theorem). The 
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verdes of T,.. were stored in san array. When the above array 


AV 


had accumulated 10 members, an approximatzon to the time 
derivation of Tay (10) - Tay (1) (divided by#unit time). The 
inclusion of the derivative feedback will give the system an 
anticipatory response. After the subtraction the array was 
updated so that a new derivative was available every 0.1 


Se@end. Due to the fact that Tay was changing at a rate of 


less than 0.5 degrees per second at distrubances below 60 


percent, the approximation of the derivative of T was con- 


AV 


Smgeredeaccurate. The sum of Tay and Te were subtracted from 


Tay references (505°F). If the error were more than 3°F an 


Serer Signal of + 1 volt was sent to the amalog cemputer. 









WORM GEAR 
ASSEMBLY 


K 
mn 
S(St_ +1) 







Parg. “(D-1). BLOCK DIAGRAM OP "Ger TEerr oO» 
POSITIONING SYSTEM. 


mWe + 1 volt error signal was placed into a standard 


K 
m 


electrical motor with transfer function S(St_ + 1) aS 


indicated in Figure (D-1) in the actual plant this motor would 
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drive at constant speed through a worm gear assembly and 
remove or insert the control rods at a constant speed. The 
Bee Seeenbly represents only a linear gaim and can be com- 
bined with the motor gain to give a new transfer function 


K 


aS The control rod effectiveness (reactivity per 
S(St_ + ‘ly 


inch) is a non-linear function of position and is different 
for individual reactors. It is assumed here that the rods 
are operating in a linear mode where 6K, = .l1 $/sec. Under 


the assumption that eae -l1 sec. the final transfer function 


1 : é 
becomes S(S + 10) ° Simulation of the drive motor yields the 


transient responses of Figures (D-2) and (D-3). From the 
Simulation of the control rod positioning device the rise 
time is 0.25 seconds and the steady state 6K, is 10.15¢/sec. 
The rise time could be shortened and the steady state ie 
remain constant if the forward gain were increased and a 
tachometer feedback loop used. However this would allow the 


possibility of an excessive 6K, 1f£ the tachometer loop were to 


R 
fail. Since the rise time of the uncompensated system was 
guite satisfactory for this application the tachometer feedback 
loop was not employed. 


Transient response of the system with external controller 


is shown in Figure (D-4) and (D-4a). 
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Pures: (D-4). TIME RESPONSE OF POWER PLANT WITH 
EXTERNAL CONTROLLER. 
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Fig.(D-4a). TIME RESPONSE OF POWER PLANT WITH 
: EXTERNAL CONTROLGER. 
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